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DENSITY ESTIMATES FOR SDES DRIVEN BY TEMPERED STABLE 

PROCESSES 

L. Huang * 1 


Abstract. We study a class of stochastic differential equations driven by a possibly tempered Levy 
process, under mild conditions on the coefficients. We prove the well-posedness of the associated 
martingale problem as well as the existence of the density of the solution. Two sided heat kernel 
estimates are given as well. Our approach is based on the Parametrix series expansion 

1991 Mathematics Subject Classification. 60H15,60H30,47G20. 

February 1, 2016. 


1. Introduction 

This paper is devoted to the study of Stochastic Differential Equations (SDEs), driven by a class of possibly 
tempered Levy processes. Specifically, we show the existence of the density, as well as some associated estimates, 
under mild assumptions on the coefficients. Weak uniqueness is also derived as a by-product of our approach. 
More precisely, we study equations with the dynamics: 

X t =x+ f F(u,X u )du + ( <j(u, X u -)dZ u , (1.1) 

Jo Jo 

where F : K + x is Lipschitz continuous, a : R + x ® K. d is measurable bounded, Holder 

continuous in space and elliptic, and (Z t ) t >q is a symmetric Levy process. We will denote by v its Levy measure 
and assume that it satisfies what we call a tempered stable domination: 

HA) < [ [ + °° 1 A (s8)^dsp(dd), VA 6 B(R d ), (1-2) 

Js*-1 Jo s + 

where q is a non increasing function, and p is a bounded measure on the sphere S d ~ 1 . Also £>(R d ) denotes 
the Borelians of R d . This is a relatively large class of Levy processes, that contains in particular the stable 
processes. 

In order to give density estimates on the solution of (ED, it is first necessary to obtain density estimates for 
the driving process. Those estimates are clear when (Z t )t>o is a Brownian motion. However, the Levy case is 
much more complicated due to the huge diversity in the class of Levy processes. Let us mention the papers of 
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Bogdan and Sztonyk [BS07I and Kaleta and Sztonyk |KS13] for density bounds concerning relatively general 
Levy processes. In the case of the symmetric stable processes, the Levy measure writes: 

r+°° r fja 

W1 e B(K d ), V (A) = / / i {seeA} C aA -—n(dO), (1.3) 

Jo Js d ~ 1 s 

for some a £ (0, 2). In the above, C a>( i is a positive constant that only depends on d and a (see Sato |Sat05j for 
its exact value), and S ' d_1 stands for the unit sphere of M d . Also, p is a symmetric finite measure on the sphere 
called the spectral measure. When the spectral measure satisfies the non-degeneracy condition: 


3C > 1, s.t. C - 1 |p| Q < f |<p,f>|“ju(d£) < C\p\ a , (1.4) 

Js d -1 

the driving process Z t has a density with respect to the Lebesgue measure. In the recent work of Watan- 
abe jWat07j . the author studied asymptotics for the density of a general stable process, and highlighted the 
importance of the spectral measure on the decay of the densities. Specifically, let us denote by pz(t, ■) the 
density of Z t , and assume that there exists 7 > 0 such that 

/i (B(0, r) n S^" 1 ) < CY 7 " 1 , V0 e S d ~ 1 , Vr < 1/2, C > 1. (1.5) 


Observe that in the case where the spectral measure has a bounded density with respect to the Lebesgue measure 
on S d ~ l , this condition is satisfied with 7 = d. For a general 7 £ [l,d] such that (11.51) holds, we have for all 


Pz(t, x) < C 


t ~ d /a 



( 1 . 6 ) 


Moreover, a similar lower bound is given for the points x £ such that two sided estimate hold in 
for 9 = x/\x\ (up to a modification of the threshold r). We refer to Theorem 1.1 in Watanabe |Wat07j for 
a thorough discussion. We would like to point out the difference between assumptions (11.51) and (11.41) . The 
assumption ED alone is enough to show the existence of the density of the driving stable process. However, it 
turns out that this sole assumption is not enough to get density estimates. Instead, we need to know what we 
refer to (with a slight abuse of language) as the ’’concentration properties” of the spectral measure to deduce 
density bounds. This concentration, reflected by the index 7 in (ED, directly impacts the decay of the density, 
as shown in the bound (11.61) . Observe however that if the concentration index 7 is too small with respect to the 
dimension, namely, a + 7 < d, the upper bound (11.61) is not homogeneous to a density, since its integral (over 
U. d ) is not defined. We refer to the work of Watanabe [Wat07 ; for a detailed presentation of these aspects. 

A generalization of this result to the case where the Levy measure does not factorize as in (11.31) . but only 
satisfies the domination m has been obtained by Sztonyk jSztlO) . Two sided estimates of the form (11.61) are 
derived, up to additional multiplicative terms involving the temperation q in ( 11 . 21 ) . with the same restrictions 
for the lower bound. 

The temperation q can be seen as a way to impose finiteness of the moments of Z (see Theorem 25.3 in 
Sato |Sat05| 1. and intuitively, the integrability properties of (Z t ) t >0 should transfer to (A t ) t > 0 . However, giving 
a density estimate on the driving process and passing it to the density of the solution of the SDE is not always 
possible. 

In the Brownian setting, if aa* is uniformly elliptic, bounded and Holder continuous, and F is Borel bounded, 
it is known that two sided Gaussian estimates hold for the density of the SDE (11.11) . This can be derived from 
the works of Friedman |Fri64] . We also mention the approach of Sheu |She91| , that also gives estimates on the 
logarithmic gradient of the density. In the stable non degenerate case, i.e. when the coefficients F, a are as 
above, and ^(d£) has a smooth strictly positive density with respect to the Lebesgue measure on the sphere, it 
can be derived from Kolokoltsov |Kol00| . that the density p(t,s,x,y) of (11.11) exists and satisfies the following 
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two sided estimates. Fix T > 0, there exists C > 1 depending on T, the coefficients and on the non degeneracy 
conditions, such that for all i,j/£ R d , 0 < t < s < T: 


c- 1 - 


- t\~ d / c 


(■ s~t ) 


(i i \ x ~v\ \ 

[ 1 + 


a+d 


<p(t,s,x,y ) < C 


- t\~ d / c 


( s~t ) 


(i i \ x ~v\ \ 

[ 1 + 


a.-\-d ‘ 


(1.7) 


This estimate is obtained using a continuity method: the parametrix technique. This approach is well suited 
to obtain density estimates for the solution of an SDE under mild assumptions on the coefficients, provided 
that good estimates are available on the driving process and on the so-called parametrix kernel. 

We refer to estimates of the form (O) as Aronson estimates: two sided bounds that reflect the nature of 
the noise of the system. In the Gaussian setting, the density of the solution has a Gaussian behavior, and in 
the stable case, the density of the solution has two sided bounds homogeneous to those of the driving stable 
process. This work aims at proving Aronson estimates when the driving process is a Levy process satisfying a 
tempered domination (in the sense of (11.21)1 . 

This research domain has a long history in the literature and is intensely developing. We can refer to 
Kochubei ' Koc89] or the book of Eidelman, Ivasyshen and Kochubei EIK04] in the parabolic setting. In the 
latter, the parametrix technique is presented as well. Other techniques to obtain density estimates comes from 
the Dirichlet froms techniques. See e.g. the works of Chen and Zhang |ZZQ13| . Our work differs from their 
in the fact that we can handle very general spectral measure with our approach. We also mention the recent 
works of Knopova and Kulik |KK14j that uses the parametrix technique for stable driven SDEs, when the jump 
part (corresponding to our coefficient a) is real valued. This assumption greatly simplifies the proof of the 
upper bound for the parametrix kernel. Also, they consider rotationally invariant stable process, which allows 
them to handle the case of a drift when a < 1. Some recent development around the parametrix in Bally and 
Kohatsu-higa |BKH14j in the Brownian setting give a probabilistic interpretation of the parametrix technique. 

Finally, we mention that existence of the density can be investigated via Malliavin calculus. In the Brownian 
setting, we refer to the works of Kusuoka and Stroock [KS841IKS851IKS87i . as well as Norris |Nor 86 l . The 
jump case is more difficult, and is treated by various authors. Let us mention Bichteler, Gravereaux and 
Jacod [ BGJ87] . and Picard |Pic96| . However, this technique requires regularity on the coefficients. In our 
approach, the convergence of the parametrix series will give us the existence of the density as well as weak 
uniqueness, under relatively mild assumptions on the coefficients. 

We will denote by [H] the following set of assumptions. These hypotheses ensure the existence of the density, 
and are those required by Sztonyk |SztlO| in order to have a two sided estimate for the driving process Z. 

[H-l] (Z t ) t >o is a symmetric Levy process without Gaussian part. We denote by v its Levy measure. There is 
a non increasing function q : R + —> R + , fj, a bounded measure on S d_1 , and a £ (0, 2), 7 £ [1, d] such that: 


HA) < 



1 A{sO)j^dsp(d9). 


( 1 . 8 ) 


We assume one of the following: 

[H-l a] p has a density with respect to the Lebesgue measure on the sphere. 

[H-lb] there exists 7 £ [l,d] such that p (B(9,r) D S d 1 ) < Cr T 1 , with 7 + a > d, and for all s > 0, there 
exists C > 0 such that: 

q{s) < Cq{2s) (1.9) 

[H-2] Denoting by tpz the Levy-Kintchine exponent of (Z t )t> 0 , there is K > 0 such that : 

E = e t¥,z(c) < e" Kt|cl “, |C| > 1. (1.10) 

[H-3] F : R+ x R d —> is Lipschitz continuous in space or measurable and bounded and a : R + x 

is bounded and ^-Holder continuous in space rj £ (0,1). 
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Also, when a < 1 , we assume F = 0. 

[H-4] <7 is uniformly elliptic. There exists n > 1 for all x,£ G K d , such that: 

<k\Z\ 2 - ( 1 . 11 ) 

[H-5] For all A G B. Borelian, we define the measure: 

v t {x, A) = v({z G a(t,x)z G A}). (1-12) 

We assume these measures to be Holder continuous with respect to the first parameter, that is, for all A G £>(R d ), 

\v t {x,A)-v t {x',A)\<C5A\x-x'\ ri{aM) f l + 1 A (s0)-^-dsp(d0). 

J s d_1 Jo s + “ 

[H-LB] There is ct non increasing function q : M-|_ —Y M_|_ and. Aiow Cl such that for all x G Ai 0 w> 

u(£(x,r)) > Ur 7 , Vr > 0 (1.13) 

Also, for all (t,x) G M+ x R d , we assume that a(t,x)Ai ow c Ai ow . 

In the rest of this paper, we will assume that [H-l] to [H-5] is in force. Also, we say that [H] holds when 
[H-l] to [H-5] hold. Note that assumption [H-2] is crucial in order to get the existence of the density. We 
point out that [H-LB] is needed for the lower bound, and that the upper bound holds independently. 

Recall that for a Markov process ( X t )t>o , the infinitesimal generator (or simply the generator) L t (x, V x ) is 
defined as the limit: 


At (a 


V x )f(x) = lim 

S^,t 


E[f(X s )\X t =x)-f(x) 


s — t 


/eC o 1 


We use a terminology from the PDE theory to emphasise the variable on which the operator acts. This speci¬ 
fication will be useful later, as we will apply these operators to transition densities, with two space arguments. 

Also, a probability measure P on Q = 2?(R + x K. d ,K.) the space of cadlag functions is a solution to the 
martingale problem associated with L t (x, X x ) if and only if for every x G for all / G Cq (M+ x K. d ,R) 
(twice continuously differentiable functions with compact support), denoting by (X t )t> o the canonical process, 
we have: 

P(X t = x) = 1 and f(s,X s ) — J^ (d u + L u (x, V x ))/(u, X u )du is a P-martingale. 

Under [H], we are able to prove the following. 

Theorem 1.1 (Weak Uniqueness). Assume [H] holds. The martingale problem associated with the generator 
L t {x,X x ) of the equation Ob 


L t (x, V x )<p(x) = ( F(t , x),V x ip(x)) + [ <p(x + aft, x)z) - ip{x) - hi v (dz), 

Jvl* 1 + \z\ 

admits a unique solution. Thus, weak uniqueness holds for EH). 

Let us mention that the weak uniqueness has been derived in similar setting by various authors. See for 
instance |Str75] or |Kom08i . We propose here an approach based on a parametrix technique, which derives from 
the works of Bass and Perkins [BP09j . 

Also, we have the following density estimate: 
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Theorem 1.2 (Density Estimates). Under [H], the unique weak solution of (ED has for every 0 < t < s, 
an absolutely continuous transition probability. Precisely, for all 0 < t < s, and x, y G R rf , 

P(X S G dy\X t = x) = p(t, s , x, y)dy. (1-14) 

Assume that the function Q defined below is decreasing, and fix a deterministic time horizon T > 0. There 
exists C\ > 1 depending on T and the parameters in [H], such that the following density estimates holds: 

VO < t < s < T, \/(x, y) G R d , 


p{t,s,x,y) < Ci 


(. s-t)~ d / a 



\y-e 3 A x )\ \ a+1 


Q{\y - 


(1.15) 


where: 

• when the drift F is bounded, 0 is the identity map: 9 s ^{x) = x, and 

— under [H-la], 7 = d and for all s > 0, Q(s) = q{s), 

— under [H-lb], for all s > 0, Q(s) = min(l, s 7_1 )g(s), 

• when the drift F is Lipschitz continuous, 9 s ,t[x) denotes the solution to the ordinary differential equation: 

■j-9 s ,t{x) = F(0 s j(x)), 9 t ' t (x) =x, VO < t < s < T, 


and 

— under [H-la], 7 = d and for all s > 0, Q(s) = min(l, s)q(s), 

— under [H-la], for all s > 0, Q(s) = min(l, s, s' 1 ~ 1 )q(s). 
Moreover, assume [H-LB] holds, and that for some C > 0, 


B(e t , T {v)-x,C(T-t) 1/a ) C A low 


there exists C 2 > 1 such that 


Co 


- t)~ d /° 


( s-t) 


1 + 


\y-0a,t(x)\ 

(s-‘) 1/a 


o+ll 


q(\y- 0 s ,t(x)\) <p(t,s,x,y). 


(1.16) 


(1.17) 


With a small abuse of language, we refer to p(t, s , x, y) as the density with respect to the Lebesgue measure. 
We point out that in [H-2], when a G R, or when p is equivalent to the Lebesgue measure on 5 d_1 this is 
actually a consequence of the Holder continuity of a, and the domination [H-l], However, in a very general 
setting, this is not a consequence of the Holder continuity of < 7 , even in the case of the stable process. Consider 
for instance Z t = (Z[,Z[), where Z\,i = 1,2 are two independent one dimensional stable processes. In that 
case, p = 6 ei + S e2 , where (ei,e 2 ) forms the canonic base of R 2 . Thus, the support of Levy measure of Z is 
exactly the reunion ejR U e 2 R. If we take aft, x) to be a rotation. There exists A G B(R . d ) and (t, x) G M+ x R d 
such that v(A) = 0, but v t {x, A) ^ 0, consequently [H-5] fails in this situation. 

The condition (11.161) appearing for the lower bound comes from the possibly unbounded feature of the 
deterministic flow associated with ED- Indeed, it states that if a neighborhood at the characteristic time 
scale of the flow stays in the sets of non degeneracy for v, then the lower bound holds. Also, in [H-LB], 
the assumption that a(t,x) G Ai ow is a form of compatibility condition. We know from Watanabe |Wat07| 
that the spectral measure plays a key role in the density estimate. This assumption ensures that the diffusion 
coefficient does not alter too much the spectral measure, so that the density estimates still holds. Besides, as 
the lower bound only holds inside the set Ai ow , we obtained Aronson estimates only in those specific cases where 
Ai ow = R d . Let us mention that [H] is satisfied in its entirety when considering a rotationally invariant stable 
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process, or a relativistic stable process. For the latter, the Fourier exponent writes (|£| 2 + 1 )“ ,/2 — 1. In that 
case, we have 7 = d, Ai ow = R d . See the remark after Theorem 2 in Sztonyk |SztlOI . 

We point out that the upper bound (11.151) is not the one established in B . rather a ’’degraded version” 
of it. Formally, it means that some integrability is required in order to perform our techniques (precisely to 
correct a bad concentration index on the parametrix kernel). Nevertheless, our approach allows to recover the 
existence and estimates on the density of m when (Z t )t >0 is a rotationally invariant stable process. See the 
proof of Proposition 13.31 and Remark 13.31 
Also, when [H-la] holds, the condition 

n(B{ 6 ,r)nS d ~ 1 ) < Cr'f~ 1 , 

actually holds with 7 = d. Besides, the function Q appearing in the upper bound (11.151) is decreasing typically 
when considering tempering function of the form: 


where m is large enough. Let us mention that when the spectral measure p has a density with respect to the 
Lebesgue measure on the sphere and the drift F is bounded, we recover the results in Kolokoltsov |Kol00j . 
without regularity assumptions. 

Remark 1.1 (On the constants). We will often use the capital letter C to denote a strictly positive constant 
that can depend on T and the set of assumptions [H] and whose value of C may change from line to line. 
Similarly, in the temperation, we will often write <7(|x|) where we actually mean g(C'|a;|). Finally, we will use 
the symbol >: to denote the equivalence: 

/ x g O 3C > 1 , C~ 1 f(x) < g(x) < Cf(x). 

The rest of this paper is organized as follows. In Section [2j we set up formally the parametrix technique, 
and give the estimates permitting the convergence of the parametrix series. Section [3] is a technical section and 
is divided in five subsections. First, in Subsection 0 we prove estimates on the frozen density. In Subsection 
we investigate the parametrix kernel and its smoothing properties. In Subsection 13.31 we tackle the well- 
posedness of the martingale problem, using estimates provided by the two previous subsections. Next, in 
Subsection 13.41 we prove the estimates giving the convergence of the parametrix series. Finally, in Subsection 
13.51 we investigate the lower bound (11.171) . 

2 . The Parametrix Setting 

We present here a continuity technique known as the Parametrix. Our approach is close to the one of Me 
Kean and Singer |MKS67) . The strategy is to approximate the solution of (11.11) by the solution of a simpler 
equation and control the distance in some sense between the two processes. First of all, let us define the proxy 
we will use. 

Let y £ R d be an arbitrary point. When F is only measurable and bounded, we take as frozen process 
X t = x + f n <t(u, y)Z u . 

Assume F to be Lipschitz continuous. Let 9 t s be the flow associated with the deterministic differential 
equation: 

-tAsOe) = F(t , 9 t s (x)), d s s (x) = x, 0 <t,s<T. 
at 

We will often refer to 0 t s (y) as the transport of y by the deterministic part of (11.11) . Recall T is the deterministic 
time horizon. Fix y £ R d , a terminal point and t £ [0, T], and x £ R rf initial time and position, we define the 
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frozen process (^J’ y )se[t,T] as: 

F(u,8 u , T (y))du + J a(u,O u ,T{y))dZ u . (2-18) 

We point out that the transport of the terminal point in the drift part comes for the unbounded character 
of the drift coefficient. Also, in diffusion coefficient a, the presence of the transport ensures the compatibility 
between the estimates on the frozen process and the parametrix kernel (see Propositions 13.11 and 13.31) . This 
technique of freezing along the curves of the deterministic system has first been introduced in Konakov Menozzi 
and Molchanov [KMM10] . and has been formalised in Delarue and Menozzi [DMIOl . In those works, the presence 
of the flow has another purpose, namely to account for the degenerate feature of their equations. However, as 
a by-product, it allows to handle the unbounded feature of the drift in a non degenerate setting. 

We could restrict ourselves to functions F that are Holder continuous. In this case, existence of the flow 9 is 
given by the Cauchy Peano theorem. However, the lack of uniqueness poses the problem of the definition of 0 tiS , 
so we decided to assume Lipschitz continuity instead. Anyhow, in the case where F is Holder continuous, we 
expect some kind of regularization by the noise, as we recover weak uniqueness, see e.g. Bafico and Baldi [BB82I . 
or Delarue and Flandoli IDF14I for recent developments. 

It is clear from the definition of (Xj’ V )sG[t,T] and assumptions [H-l] (domination of the Levy measure) and 
[H-2 ] (non degeneracy of the spectral measure) that X T,y has a density with respect to the Lebesgue measure. 
We denote the latter: 

P(Xj ,v £ dz\Xj’ v = x) = p T,v (t , s, x, z)dz, s £ (f, T}. 

To get an explicit representation for it, we proceed by a Fourier inversion. Observe first that the Fourier 
transform of Xj' y actually writes: 

= e i<P,x+J t a F(u,e u , T (y))du) exp lpz ( CT (u i e UtT {y))*p) 

where we denoted by a(u,0 Ut T(y))* the transpose of a(u,9 Ut T{y)), and tpz is the Levy-Kintchine exponent of 
Z. Due to assumptions [H-l] and [H-2], this exponent is integrablc and the frozen density actually writes: 



p T - y {t,s,x,z) = 


(2tt)‘ 


■ / dpe~ i< ' p,z ~ x ~^ F (uJv.,T(y))du) 

Js. d 


x exp ( / du e 


(p,‘ 


_ ! _ , (p,g(M u ,r(y))Q . 


l+l £| 5 


<d0 


(2.19) 


We will often denote p(t,T,x,y) = p T,v (t,T : x,y ), namely, we omit the superscript T, y when the freezing 
parameters and the points where the density is considered are the same. Observe that in this case, we have 


y- J F(u,0 UyT (y))du - x = 8 t , T (y) - x. 


( 2 . 20 ) 


Proposition 2.1. Assume that the semigroup (P s ,t)o<s,t<T generated by the solution (A s ) sg r t r] to (11.11) is 
Feller. We have the following formal representation. For all t > 0, (x,y) £ (R d ) 2 and any bounded measurable 
f : R d -a- R: 

P T ,tf(x)=E[f(X T )\X t =x}= f f(p«flW)(t,T, I ,y) f(y)dy, (2.21) 

\r=0 J 

where Ft is the parametrix kernel: 

VO < t < s < T, (x, y) £ (R d ) 2 , H(t, T, x, y ) := (L t (x, X7 x ) - L t (9 t , T (y),X x ))p T ’ y (t, T, x, y). 


( 2 . 22 ) 
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The notation ® stands for the time space convolution: 

f ®g(t,T,x,y)= l du dzf(t,u,x,z)g(u,T,z,y). 

Jt J R d 

Besides, p® = p and Vr £ N, H^ ( t , T, x, y) = H <g> H(t, T, x, y). 

Furthermore, when the above representation can be justified, it yields the existence as well as a representation 
for the density of the initial process. Namely ¥[Xt £ dy\X t = x] = p(t,T,x,y)dy where : 

+ OO 

Vi>0, (a :,y) e (R d ) 2 , p(t,T,x,y) =^2(p® H (r) )(t,T,x,y). (2.23) 

r —0 


Proof. Let us first emphasize that the density p T ' v ( t , s, x, z ) of Xl ,x ' T ’ y at point 2 solves the Kolmogorov backward 
equation: 


dpTv 

dt 


(' t,s,x,z ) 


-Lt{0t, T {y),Vx)p T,v (t,s,x,z), for all t < s , (x,z) 


p nd 


x R nd , lim t\sP T ’ v {t, S, •, z) 


«*(■)■ 

(2.24) 


Here, the notations V x highlights the fact that L t (9 t ,T{y), V^) acts on the variable x. Let us now introduce the 
family of operators ( Pt,s)o<t<s ■ For 0 < t < s and any bounded measurable function / : M. nd —>■ K: 


Pt,sf{x) 


■= [ P{t,T,x,y)f{y)dy 

JR nd 


p T ’ v (t, T, x, y)f (y)dy. 


(2.25) 


Observe that the family ( Pt, s )o<t<s is not a two-parameter semigroup, because of the integration with respect 
to the freezing parameter y. Anyhow, we can still establish, see Lemma 13.21 that for a continuous /: 


lim P s ,tf(x) = f{x). 


(2.26) 


This convergence is not a direct consequence of the bounded convergence theorem since the freezing parameter 
is also the integration variable. It rather follows from the estimates on the Frozen density p T ’ v (t,T,x,y). 

The boundary condition (12.2611 and the Feller property yield: 

(■ Pt,T - Pt,r)f(x) = jt du-^-^P tlU (P u ,Tf(x))^. 

Computing the derivative under the integral leads to: 


(P t ,T ~ P t<T )f(x) = J du^d u P t , u {Pu,Tf(x))+Pt,u{du{Pu,Tf{x)))j. 

Using the Kolmogorov equation (12.24[l and the Chapman-Kolmogorov relation d u Pt,uP(x) = Pt, u {L u (x , V x )(p(x)), Vy> G 
l) we get: 


r 2 


(Pt,r ~ Pt,r)f{x) = Jt duP t>u (L u (x,X x )P UtT f'){x) - P t , u \^J J{y)L u {9 u ^{y)y x )p(u,T,-,y)dyj (x). 

Define now the operator: 

'Hu,t(<p)(z) ■= [ T(y)( L u(x, V x ) - L u (9 UiT (y),X x ))p(u,T, z,y)dy = [ ip(y)H(u,T, z,y)dy. (2.27) 
J R" d ' ' Jw nd 
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We can thus rewrite: 

Pt,rf(x ) = P t , T /( x) + J P t ,u('Hu,T(f)){x)du. 

The idea is now to reproduce this procedure for P ti „ applied to TL u ,T(f)- This recursively yields the formal 
representation: 


_ _ pT pU i r-Ui —i 

Pt,Tf(x) = Pt,Tf(x) + ^2 du i / du 2 ... du r Pt,u r ('H-u r ,u r - 1 o ■ ■ ■ oH Ul , T )(f)(x). 

r>1 Jt Jt Jt 

Equation (12.211) then formally follows from the following identification. For all r G N*: 

pT pU\ pu r — i p 

du i / du 2 ... j du r P t ,u r {'Hu r ,u r - 1 °---°'Hu 1 ,T)(f)(x)du= j f{y)p®H ( ' r \t,T,x,y)dy. 

ftJt 


We can proceed by immediate induction: 

/ T pu i ru r — i 

dui J du 2 . .. J du r Pt, Ur ( / Hu r .,u r - 1 °---°'Hu 1 ,T){f)(x)du 

nT pU\ pu r - 1 p 

= / du! du 2 ... du r dz'Hu r ^ur-iO--- 0 'Pu 1 ,T{f){z)p{t,Ur,x,z) ( 2 . 28 ) 

Jt Jt Jt J R nd 

pT PU 1 ru r -1 P P 

/ du! du 2 ... du r dz dyH Ur _ ljUr -2 ° '' ’ ° 'H Ul ,T{f){y)H(u r , u r _i, 2:, y)p(t, u r , x, z) 

Jt Jt Jt J R nd J R nd 

pT pU\ pu r -2 P 

/ du1 l du 2 ... dur-i dyH Ur _ ltUr - 2 0 • • • 0 'H Ul ,T{f){y)p a ® H(t,u r -i,x,y). 

It Jt Jt Jm nd 


(t2~2Tl 


Thus, we can iterate the procedure from (12.281) with p® H instead of p. 


□ 


The existence of the density for the solution of ed will follow from the convergence of the parametrix series. 
In the following, we will denote 


p(t,T,x,y) 


(: T-t)- d / a 



|y-6> Tif (x)[\ a+7 

( T-t )V“ j 


q(\y - d T ,t(x)\). 


(2.29) 


This is the upper bound on the Frozen density under [H] derived by Sztonyk ISztlOj . adapted to our possible 
unbounded drift case. We prove that this upper bound holds for the frozen density in Section [3] 

The following lemma proves the convergence of the series (12.281) . 

Lemma 2.2 (Control of the iterated kernels). There exist ( pn\ = QrZl T ) > 0, w G (0, l] s.t. for all 

t€[0,T], (x,y) G (R d ) 2 : 


< 


< 


322 )( (T _ 


((T ~ t)“p(t, T, x, y) + p(t, T, x, y)), 


qm (T 


(T — t) u p(t,T,x,y), 


\p®H(t,T,x,y)\ 
\p <g> H(t,T, x, y)\ 


UJ 


(2.30) 

(2.31) 
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where we denoted p(t , T, x, y) = S A \x — d t ^T{y)\ v< ' aM ' ) p{t, T, x, y). Now for all k > 1, 

I V ® T , x,y) | < (q^) 2fc ( ^r (( T - *)*"£(*> ^ ^ v) + (p + p)(t, T, x, y)) , (2.32) 

|p ® (2fe+1) (*, T, X, y ) I < (C^]) 2fc+1 (V - *) (fc+1)w P + (T - tf (p + p)+ p) (i, T, x, y), (2.33) 

where we set u>' = uilT-t <l + kixlT-t>i > 0. 

The above controls allow to derive under the sole assumption [H] the convergence of the parametrix series 
(thus, existence of the density for the solution of (11.11) 1. and the upper bound (11.151) for the sum of the parametrix 
series (12.231) . 


Proof. Let us denote by 


H{t, T, x, y) 


Sa \x - d t , T (y)\ ri{aA1) 

T-t 


p(t, T, x, y). 


We prove in Section EQ1 that H is actually the upper bound for the kernel H. The proof of Lemma [321 relies 
on the two following estimates: 


p(t,u,x, z)H(u,T, z,y)dz < Cq((T — uf 0 1 + (u — t) u i + 


„_i , (5 A|x-0 t ,r(y)r (QA1) \_ 


T-t 


^Jp(t,T,x,y) 


j*p(t,u,x,z)H(u,T,z,y)dz <C 0 ({T - u) u 1 + (u-t) u l ^jp(t,T,x,y). 


(2.34) 

(2.35) 


Integrating these estimates in time yields (12.301) and (12.311) . setting = 2Cq. We postpone the proof of 
these important estimates in Section [3j as it is technical and relies on sharp estimates on the frozen density and 
on the parametrix kernel (see Lemmas 13.61 and 13.71) . Assuming estimates (12.341) and (12.35|) . we prove estimates 
(12.321) and (12.331) by induction. 

Initialization: 

Since (T — t)“ (p + p) > 0, we clearly have: 


I P ®H(t, T, x,y )| < Cfosifc T - t) u p + p+(T -t) u (p + p)') (t, T, x, y). 


We now turn to the estimate for | p 0 H^\. Starting from (12.301) . we use equations (12.341) and (12.351) . we have 
to write: 
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/^jt_- r 1 /» 

\p <g> H^ 2 \t, T, x, y)| < -&21 / du(u — t) u / p(t,u,x,z) 

U Jt JWL d 


5 A\x - 0 u , T (.y)\ v{aM] 


T-u 


p(u,T,z,y)dz 


< 


,<te \ T, f u A |x — 0 uT (y)\ r,( - aA1 ^ m 

+->— 1 / du / p(t,u,x,z J-—- p(u,T,z, y)dz 

u Jt Jm d T ~u 

jT du(u - tr ((T - ur- 1 + (u- t ) w_1 + 5 ^ X ~ 9 ^ T }f lV{aM) )p(t, T, x , y) 
+ ^C 0 ^ dw(T - u )“ _1 + (u - t)^ 1 ^ p(i, T, x , y) 


< 


< 


CO 

9 m 

Ld 


Co 


(-B(w + 1, w) + (T - t) 2u p(t , T, x , y) + ^ p(t, T, x, y) + 2 ^ p(t , T, x, y) 

(T - t) w ((T - i)“p + (p + y)) (i, T, x, y). 


To get to the last estimate, observe that the biggest contribution amongst the constant multipliers is We 
also recall that Ofc%[ = 2C 0 . 

Induction: 

Suppose that the estimate for 2k holds. Let us prove the estimate for 2k + 1. Recalling that p and H are the 
upper bounds for p and H respectively, we write: 


\p (g> ij( 2fc+1 ) (t, T, x, y)\ 

< — ( f du[u — t) 2ku+l f p(t,u,x,z)H(u,T,z,y)dz + f du(u — t) ku f (p + p)(t,u,x, z)H(u,T, z,y)dz 

jR d Jt J«. d 

£ du{u - t) 2ku C 0 ((T - ur - 1 + (« - tr- 1 + 6 *\ x -^vW laA1) y )p(t, T, x, y) 

S A \x-e t ,T{y)\ v{aM) ' 


< 


< 


k\co 2k 

(teD 2fc fT 

k\co 2k +' Jt 

+ ^S£- [ T du{U - t)kU1 °° ( 2((T “ Mr_1 + (U - * )W_1 ) + dAlX ~^ )r ^ ; )pr T, x, y) 

(qo) 2fc 

k\co 2k 

i (qfczr 

k\co 2k 

(qo ) 2fc+1 


( B<2t “ + '• “> + (STIv) (T - t)<2 ‘ +1, "« i ’ r ' *■ »> + T - *■»> 

2^B(ku) + 1, w) + ^ fc+ 1 - L ^ ) ( T - t) (fe+1)w p(t, T, x, y) + ^ ^ y(t, T, z, y) 
(fc + i^fc+i < T - *)*“ (( T - *) (fe+1) ^ + (T - tf (P + p)+ P ) (t, T, x, y ). 


Co 

Co 


To get the last inequality, observe that from all the constant multipliers, the biggest one is 
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Assume now that the estimate for 2k + 1 holds. We prove the estimate for 2k + 2 . 
\p®H( 2k+2 \t,T,x,y)\ 

J\2fc+1 / r T 


< 


(tel ) 2 


k\u 2k+1 

rT 


f du(u — t)( 2fc+1 )“ f p(t,u,x,z)H(u,T,z,y)dz 

Jt J R d 


+ / du(u — t ) 




{p + p) (t, u, x , z)H(u, T , z, y)cb 


+ / du(u — t) ku / p(t,u,x,z)H(u,T,z,y)dz 

Jt Js. d 

.\2fc+l /-T 


< 


(tei ) 2 


£ du(u - t)<- 2k+ V“Co ((T - u )"" 1 + (« - t )"" 1 + T) x, y) 

£ du(u - t) k “+“'C 0 (2 ((T - U )"" 1 + (u - t)"" 1 ) + * A I g-gt,r(y)r (aA1) )- (t] T) y) 

J" T du(u - t) k “C 0 ((T - u )"- 1 + (« - i)— 1 )#*, T, x, y) 

c “ te-+ i >-+u-)+^ *■»> +r - *■ »>) 


+ 


< 


fc!w 2fe +! 

(te 2fc+1 

fc!a; 2 fe + 1 

1 (tef fc+1 ' r 

fc!a;2fe+i 

(te 2fc+1 

/ c ! w 2 fe + 1 

, (teP 2fc+1 


Cn 


(2k + 2)u> 
2 


jfe!w 2 *+ 1 

(teD 2fc+1 


( (2B(fcu; + w' + 1, w) + TT-r^-— 7 ) ( T ~ t) (fe+ 1 ) “+ w 'y(t, T, x, y) 
\\ (fc+ljCJ + u// 


(2fc + l)a; 

_ £^fcu;+u/ 

/ccj -|- Ct/ -|- 1 


p{t,T,x,y) 


< 


k\uj 2fe + 1 

(tef fc+2 

(fc + l)!w 2fe+2 


Co ( B(kto + l,w) + 




(fc + l)w, 

(T _ t) (fc+i)« (( T _ t )(fc+i)«-p( t) T, x, y) + (p + p)(t, T, x, y)) . 


To obtain the last inequality, and the estimate for p <g) H 2 ( k+1 \ we factorised by (T — f)( fc + 1 ) w . Also, the 
biggest constant multiplier is prpryy- 

□ 


3 . Proof of the estimates. 

In order for the parametrix technique to be successful, we must obtain some sharp estimates on the quantities 
involved in the parametrix expansion (12.231) . This is usually done in two parts, first, we give two sided estimates 
on the density of the frozen process, as well as a similar upper bound on the parametrix kernel H , up to a time 
singularity. Then, we prove that those bounds yield a smoothing effect in time for the time space convolution 
p® H appearing in (|2.23l) . 


3.1. Estimates on the Frozen Density 

We first give the estimates on the frozen density. 

Proposition 3.1. Assume [H] is in force. There exists C > 1 s.t. for all t £ [0, T], (x,y) £ (R d ) 2 : 


p T ’ v (t , s,x, z) < C- 




(s-t) 


\z-x-ff F(u,e n , T (y))du |\“+ 7 


1 + 






q C 


-1 


z-x- F(u,0 U ' T (y))du 


(3.36) 
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Moreover, assume [H-LB] holds. Then, for all z — x — J f s F(u, d u ,T{y))du £ Ai ow , the lower bound holds: 

C l ~ 1 (s — t)~ l 

If F{u,O u<t 
(s-ty/- 

Proof. We prove these estimates in the lines of Sztonyk [SztlOj . The idea consist in splitting large jumps and 
small jumps at the characteristic time scale and exploit the Levy-Ito decomposition. Fix t,T £ R + and y £ R d . 
Observe that since the drift part in the frozen process is deterministic, it suffices to prove the estimates or: 

A s = a(u , O u , T {y))dZ u . 

We point out that A s = A s (t, T , y), where t, T , y are fixed. The Fourier transform of A s writes: 

E ( e dC, A d) = exp du J ^ e l<C ’ CT( “’ e “' T(y))c> - 1 - i((, a(u, 0 lt ,T(y))£) 1 {|£|<iV«}K<^ ■ 

Changing variables in the time integral to v £ [0,1] and setting a v = cr((s — t)v + t, 0( s -t)v+t,T(y)), we obtain: 

E(e l<c ’ As >) = exp ({s-t) ( du [ e t < c ’ av ® - 1 - i{(, 

\ JO JR d 

Now, defining vg to be the image measure of dw(df) by the application (v,£) i —> a v f, we obtain: 

E(e l<C ’ As> ) = exp (is - t) e l{c ’ v) - 1 - i{C,rj}l {M < t i/c } ^g(dr]) S j , 

which is the Fourier transform of some Levy process (<S„) U >o, with Levy measure vg, at time s — t. In other 
words, the marginals of (A s ) sg r t T i corresponds to the marginals of (S s _t)se[t,T] : 



i/a. 


{y))du\ 




q[C 


z — x — 


F(u,6 u , T {y))du^j <p T ’ v {t,s,x,z). (3.37) 


Vs G [t,T\, A s 


(law) 

— '-’s-t- 


(3.38) 


The idea is now to work with the process (<S u ) u >o, and prove that it satisfies the assumptions [H]. Specifically, 
we prove that [H-l] and [H-2] holds for vg, and that when [H-LB] holds for v, it holds as well for vg. 

Let A £ B( R d ). By definition of vg , we have: 

vs{A) = / / 1 {<r„t£A}v{d£)dv. 

Jo J R d 

From the tempered stable domination, we deduce: 

vg(A) < J J ifcr v s^A}^^-fJ'((k)dv. 

For fixed v,q, we change the variables in the integral in ds to p = s|<r„?|. Observe that from the uniform 
ellipticity of a and the doubling property of q, we have q < Cq(p). It yields: 


”s{A) < 


■ 

Jo Js d ~ 1 Jo 


. _ <i{p)dp 

pl+a 


\a v q\ a p,(d<;)dv. 


+00 
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We now define ^s(efc) to be the image measure of \a v <;\ a p{d^)dv (measure on [0,1] x S d *) 
(v,<;) i-A jfjdy- We thus obtain: 


by the application 


vs (A) < 


S d_1 Jo 


+ °° q(p)dp 

1 {pCeA> , +a Ps(dQ. 


Consequently, [H-l] holds for v$. Observe that by construction, we have that [H-2] holds for ps- Therefore, 
from Stzonyk |SztlO] , denoting ps(u, •) the density of S u , the following upper bound holds: 


Ps{u : z) < C 


l~ d / a 


(i + *) 


a +7 


3(1*1)- 


We deduce the estimate for A s , and the upper bound on p T ' y (t, s, x, z) then follows. To see that the constant C 
above does not depends on t, T, y, we have to see that ps(B(x, r ) 0 5 d_1 ) < Cr 7-1 , where C does not depends 
on f, T, y. By definition of ps, we have: 


Ps(B{x,r) nS d ') = j r)n sd-i}Wv<;\ a Kd^)dv. 

Now, observe that £ B(x,r) => ? £ B(Ca~ 1 x,Cr), where C is the uniform ellipticity constant of a(t,x). 
Consequently, 

ps(B(x,r) H S^ 1 ) < f p[B(Ccr~ 1 x, Cr) n S d ~ x )dv < Cr 1 ~ 1 , 

Jo 

uniformly in the parameters t,T and y. To get a lower bound on p T ’ y (t,s,x,z), we investigate a lower bound 
for ps(u , •). To that aim, we prove that when [H-LB] holds for v, it does for vs- Specifically, assume [H-LB] 
holds for v. By definition of ps, for all x £ R d , r > 0, we have: 


vs 



^-{\x—ct v ^\ <r}dw [d<J). 


Now, from uniform ellipticity of cr, 


{\x - cr„?| < r) D {|cr„ X X - ?| < Cr}. 


Now, since by assumptions, a v 1 x £ Ai ow for all v £ [0,1], we have: 


i/ s (B(x,r)J = f v(a~ 1 x,Cr'jdv > j Cr< dv>Cr^fi X ^ 


Cv X 


a +7 


a+7 ’ 


where to get the last inequality, we exploited the uniform ellipticity of cr. Besides, for r £ (0,1), we write 
using the ellipticity of cr: 


v s (B(0,r) c ) = [ I 1 { | CT „ ? | > r }dvv(d<;)< j [ 1{| q \>Cr}dvv(d<;) < c\. 
v 'Jo J R d Jo JWL d r 

Thus, we recovered [H-LB] for vs and the lower bound holds for ps(u,-). Thus, the one for p T,y (t,s,x,z) 
follows. 

□ 
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Remark 3.1. The idea of the proof was to identify the density of the frozen process to the density of some 
Levy process and exploit the Levy structure to derive bounds on the density. The procedure described above 
require the uniform ellipticity of a in order to prove that the assumptions [H] holds for the new Levy process 
(iS„) u >o- Intuitively, we can say that a uniform elliptic coefficient does not alter much the nature of the noise in 
the system. From the identity in law (13.3811 that holds for fixed s £ [i, T] and equation (12.201) . we deduce that: 

p T ’ v (t,T,x,y) = ps(T - t,0 t ,T(y) - x). (3.39) 

This identity will be useful when investigating the parametrix kernel H. 

Now, we state a Dirac convergence Lemma for the frozen process when the freezing parameter changes. This 
convergence will be used in the proof of the well posedness of the martingale problem. The difficulty comes 
from the fact that when integrating with respect to the freezing parameter (as it is the case in a parametrix 
procedure), the Dirac convergence does not follow from the Chapman-Kolmogorov equations. However, since 
we have good estimates on the frozen density, we manage to prove the following lemma: 

Lemma 3.2. For all bounded continuous function f : —>• R, x £ R d , 


[ f(y)p T ' v (ti T, x, y)dy — f{x) 

J R d 

that is, for all (x,y) £ R d x p T,v (t,T,x,y)dy =>■ S x (dy) weakly when T ft. 
Proof. We prove this convergence in the lines of | IHM14j . Let us write: 



(3.40) 


[ f(y)p T ' v (f,T,x, y)dy — f(x) = [ f(y)(p T ’ v (t,T,x,y) - p r ’ 0T ’ t< - x \t,T,x,y))dy 

JR d JR d v 7 

+ [ f(y)(p T ' 9T ’ t( ' x) {t,T,x,y))dy - f(x). 

J R<* V 7 

From the usual Dirac convergence in the Kolmogorov equation (12.241) . the second term tends to zero when 
T —» t. We focus on the first term. Define: 

A = f f{y)(p T ’ v (t, T ,x,y) - p T ’ eT -^ x) (t,T,x,y))dy. (3.41) 

J r<* v 7 


For a given threshold K > 0 and a certain (small) /3 > 0 to be specified, we split R. d into D\ U D 2 where: 


D\ = 


{y e !/«' - K{T ~ tr> > D * = {y e |gt,T ffa x| > K( ?-v~ 0 } • 


A direct application of Proposition 13.II yields: 


' _ f\-d/a 


p T ’ v (t, T, x, y ) < C- -^— — — —^q(\e t>T (y) ~ ar|). 

/1 , Wt,T(y)-x\\ 

V + (r-t)» 7 - J 


Now, observe that 


y-x- 


f F(u,6 u>T {9 Tjt (x)))du = y-x- f F(u,9 Uit (x))du = y - 0 T , t (x). 

Jt Jt 
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Also, from the Lipschitz property of the flow, we have 1 0t,T(y) — x\ x \y — dT,t(x)\. Consequently, we obtain: 

ptfi T Ax)(t,T,x,y)<C (r ~^ d/a 


and we have the same upper bound for the two densities in (13.411) . The idea is that on £> 2 , we use the tail 
estimate, and on D 1 , we will explicitly exploit the compatibility between the spectral measures and the Fourier 
transform in the Fourier representation of the densities. Set for i € {1, 2}: 

Ad, ’■= J f(y)(p T ’ v {t,T,x,y) -p r ’ 0T ' t{x) (t,T,x,y)')dy. 

For D 2 , we bound the two densities as we described above: 


1 + 


\0t,T(y)-x\ 

(T-t)V“ 


a+7 


?(l 9t, T {y) -x|), 


\!d 2 \ < 

C\f\oo 

/ ^ * } d/ ° - *|)dy 

A (1 + 



r +00 r d-l 

< 

C\f\oo 

/ , a+ ~Q(r)dr 

J K(T—t)~P 1 + r +7 

< 

C(T — 

^(■y+ot-d) 


Thus, for /? > 0, Id 2 —> 0. On D 1 , we will start from the inverse Fourier representation of p T,z (t,x,y ), 

2 = dT,t(x),y. Recall we denoted by ip z the Levy Khintchine exponent of Z, that is e t(pz ^ = E(e^ p,Zt ^), 
denoting cr* the transpose of tr, we have: 


p T ’*(t,T,x,y) = -±- d [ d( e -«C,y-L T exp ( F <p z ( a (u, 6 u , T (z))*()du ). 

{ 2 n) a J R d \J t / 


Consequently, we have: 


p T ’y(t,T,x,y)-p T ’ eT ^ x \t,T,x,y) 


1_ f e -*(C.y-/ t T F ( u ,8u, T (y))du-x) e f t T vz(o-(u,6 U ' T (y))*()du 

n ) d J R“ 


(27TJ“ J R d 

y~St F i u ^u,t(oc))du-x) e f^ <pz(<r{u,9 u ,t{ x )YC) du ( lQ 


(2 


1 _ f (g-id’V-f? F ( u fiu,T{y))du-x) _ e -i{C,y-Jt F {u,e u ,t{x))du-x)\ e St Vz(<?{u,9u, T (y))*O du c lf 

' ) 

I ^ f g-*(C,y-/ t T F (u,9 n , t (x))du-x) <pz(<r(u,9 U ' T (y))*O d u _ e f? <pz(cr(u,9v„t(x))*Odu\ ^ 

(27r) d J R d V ) 

= T 1 (t,T,x,y) + T 2 (t,T,x,y). 


We split accordingly: 



f{y) ( P T ' v {t , T, a;, y) - p T ' eT ’* (x) (t, T, z, y)) dy 



f(y)Ti(t,T,x,y)dy 



f(y)^ 2 (t,T,x,y)dy. 
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Note first that when a < 1, we assumed F = 0, so that the term Ti(t, T, x,y) = 0 in that case. We now treat 
this term, with a > 1. Using the mean value theorem, we write: 


ri (t, T, x, y) 



where we denoted by I the identity map of Recall that from the Lipschitz property of the flow and Gronwall’s 
Lemma, there exists C > 0 such that for all t < T, z £ R d , | (I — 0T,t)(z) | < C(T — i)(l + |z|). Thus, since 
y £ D i, we have for fj < 1 /a , 


x,y)\ < C(T-t) [ Ide-^-^dC^CiT-t) 1 -^. 

Js. d 


Integrating on D±, we obtain: 


f{y)Fi(t,T,x,y)dy 


>Dl 


< C\f\oo{T — t ) 1 


T-yt 


when l/d(l — 1/a) > /?. For T 2 , we write: 


r 2 (t, T, x, y) = 1 f d£ e - i (C,y-0 T ,t(.x)) f ^ x vz(<?(u,e u , T (y))* 0 +(i-\)'pz(<T(u, 8 u,t(x))*(;)du 

(2ir) d J Rd J 0 

x £ (Fz(cr(u, d UtT (y))*C) - t Pz{<r{O u ,t{x))*C))du. 

We know from assumption [H-2] that the Levy-Khintchine exponent is bounded by — K(T — t)|CI“> thus, we 
obtain independently of A £ (0,1): 

e / t T Ay?z(er(u,0 UiT (y))* 0 + (l-\)<pz(<r(u,6 Ui t(x))* ()du e ~K(T~t) |C|“ 

On the other hand, using the bound on the Levy-Khintchine exponent and assumption [H-5], we can rewrite 
the increment: 

ipz(cr(u, 9 u , T {y))*0 - V>z{a{u, 9 u ,t{x))*C}d 
= [ [ cos((a(u,6 u , T {y))*C,0) - cos((a(u,e u! t(x))*</,^))u(dz)d 

Jt J K d 

< K\c\ a ^ \e Uyt (x) - 6»„, T (y)r (aA1) < K{T -1) |C|°|* - 0 t ,T(y)r (aA1) . 

To summarize, we obtained: 

[ f(y)T 2 (t,T,x,y) < 1/U [ dy\T 2 (t,T,x,y)\ 

J Di JDx 

< cf dy [ (T - t)|C|“|* - 
J Di JRd 
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Changing variables, and integrating over f yields 





< 


C 

t d /° 


l/lc 


id 


C 


dy\e t , T {y) ~ *\ v{aA1) 


F-tyvM-Jo 


d r r v(aiAl)+d-l _ ^d./a+r)(lM/ot)' 


Choosing now > /3 > 0 gives that | Id 1 | —> 0, which concludes the proof. 


T\t 


□ 


3.2. The Smoothing Properties of H(t,T,x,y). 


First, we investigate an upper bound for the parametrix kernel. Then, we use it to prove a smoothing effect 
in time, in the sense that the singularity in time is relaxed after an integration in space. Recall that: 

Vt > 0 , (x, y) £ (R d ) 2 , H(t, T, x, y) := (l{x, V.) - L{G t , T (y), V x ))p T ’ y (t, T, x, y). 


Heuristically, the difference of the generators should yield a singularity (T — i) -1 , just as in the Gaussian case. 
This singularity is not integrable in time, however, we will prove that after an integration in space (which 
naturally occurs in the time-space convolution in the parametrix series), it will be. 

Proposition 3.3. Assume [H] is in force. There exists C > 0 s.t. for all t £ (0,T], (x,y) £ (R d ) 2 : 


\H(t,T,x,y)\<C((T-t) 


W 1 /" ■ 


{«>!} 


6 /\\x-6t, T {y)\ r,{aA1) 


T-t 


p{t,T,x,y) = H(t,T,x,y), 


where 

• when the drift F is bounded, 6 is the identity map: 0 t j^(y) = y, and 

— under [H-la], 7 = d and for all s > 0, Q(s) = q(s), 

— under [H-lb], for all s > 0, Q(s ) = min(l, s 7_1 )q(s), 

• when the drift F is Lipschitz continuous, d s ^{x) denotes the solution to the ordinary differential equation: 

-rOs.riy) = F{0 s , T (y)), 9 t,t(jj) = y , VO < s < T, 


and 

— under [H-la], 7 = d and for all s > 0, Q(s) = min(l, s)q(s), 

— under [H-la], for all s > 0, Q(s) = min(l, s, s 7 _ 1 )q(s). 

Thus, the upper bound on the Kernel H is the same as the upper bound on the Frozen density p T ’ v (t , T, x, y) 
up to the additional multiplier (T — t) - 1 /“l a >i + (S A \x — 0 t ,T{y )| r, (“ A1 ))(T — t) -1 , that can be seen as the 
singularity induced by the difference L(x, V x ) — L(6 t ,T(y), V x ) applied to the frozen density. The proof proceeds 
following the lines of Sztonyk |SztlO| . splitting the large jumps and the small jumps. The small jumps are dealt 
using Fourier analysis techniques, whereas the big jumps are dealt more directly. 

Proof. Recall that from (13.381) . the density of Xj’ y can be linked to the density of the Levy process (S u ) u >0 
considered at time s—t. We now exploit the Levy structure of (S u ) u >0 to obtain an upper bound on H(t, T, x, y). 
Specifically, let us introduce the Levy-Ito decomposition of (S u ) u > 0 : 

S u = M u + N u , 

where (M u ) u >0 is a martingale and (N u ) u > 0 is a compound Poisson process. We choose to place the cut-off at 
the characteristic time-scale, namely (T — t) 1 / 01 . Therefore, the Fourier transform of M u writes: 
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E (e*<C,M u >) = exp (u J^( e *(C,ri) _ l-i(C,?j))l { |,|< (T _ t) i/a } J/s(d^ ■ 

This expression is integrable and regular in the variable £ (see Section 2 in Sztonyk |SztlO| and the references 
therein). Thus, the density pm(u , •) of M u exists and is the Schwartz’s class. Thus, we can say that this term 
produces the density in the Levy-Ito decomposition. Also, denoting by vsidz) = ^-{\z\>(T-ty/ a } v s{dz), we have 
the following decomposition for the law of the compound Poisson process N u : 


Pn u ( dz ) = e 


+ OO 


k=0 


,k r,*k 


(dz) 


k\ 


Now, by independence of (M u ) u >o and (N u ) u >o and exploiting equation (13.391) . we get: 

p T ’ v (t,T,x, y) = p s (T — t, 9 t , T (y) — x) = [ p M (T - t, d t , T (y) ~ x - QP NT _ t (d£)- (3-42) 

JR d 

From the definition of the generators, the operator naturally splits into three parts, for a test function ip: 


( L t {x , V x ) - L t (Q t ' T (y), V*)) tp(x) = (V x ip(x), F(t, x) - F(t , 0 t>T (y ))) 
+ J (<p(x + z)~ v(x) - (V^(x), z)^j l { | z |<(T-t)i/°}(i / t(x, dz) - v t (d t ,T(y),dz)) 
+ J (f(x + z)~ p(x)^j l{| z |>(T-t)i/«}(^(x, dz) - v(0 t , T (y), dz))- 


Recall that we defined Pt(£,A) = v{z £ K d \cr(t,£)z £ A}. Also, observe that by symmetry of v, we changed 
the cut-off function to exhibit the intrinsic time-scale. Note that the first order term in the operator is present 
only in the case a > 1. Otherwise, we assumed that F = 0. 

Thus, a derivative along x of p(t , T, x, y) acts in fact on the density of the martingale, and we have to control 
V x pm(T — t, 0t,r(y) — x — £). Borrowing the notations of the proof of Lemma 2 in [SztlO] . we have: 

pm(T -1,0 t , T (y) -x-0 = (T- t)- d / a g T -t((T - ty^e^iy) -x-€). 

Formally, since we chose to split at the characteristic time-scale (T — f) 1//a , the density of the martingale 
presents a time space separation, and defining gr-t as above allows to have estimates independent of T — t. 
Thus, uniformly for all T — t > 0, gT-t(-) is in Schwartz’s class. Therefore, we have for all m > 1: 


I V*Pm(T - t, 0 t , T (y) ■ x £)|< {T _\y, a C m (T - t)~ d / a (l + ^ T ^ ] _ f /a ^ 

and we recovered Lemma 2 in ISztlOj . up to the singularity 1 /(T — t). Consequently, integrating this estimate 
against the law of the large jumps, we obtain the following estimate on the gradient of the density: 



S7 x p(t,T,x,y) < 


C 


T-t 


(T-t)~ d / a 



|flt,T(y)-zi \ Q+7 
(T-tp/- ) 


y\0t,r(y) 


x|) 


Thus, a derivative on the density yields a singularity in t 1 ' a which is integrable when a > 1. Specifically, 
when the drift F is bounded and a > 1, we write: 
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I<v x p(i, T, x, y),F(x) - F(0 t , T (y))}| < C f 2| J P| 00 |V sc p(* ) T, x, y)\ 

< c(t- ty'/’ {T -‘ ) -“ / ° mie,, T (v) - *D- 

( 1 + lteM^l) 

On the other hand, when F is unbounded, we have to deteriorate the tempering function, exploiting the 
Lipschitz property of F: 


\(V x p{t,T,x,y),F{x) - F{0 t , T (y))}\ < C\x — 0t, T (y)\\V x p(t, T, x, y)\ 

' - t\~ d / a 


< C(T-t)~ 1/a - ——— 


1 + 


\ s t,T(y)-x\ ^j a+7 


< C(T-t)~ 1,c 


(T-t) 

(T -t)~ d /° 


(1 I |gt,T(y)-a-Q 

V + ( T-ty/o, ) 


a+7 


\x - 0t, T (y)\q(\O t ,T(y) -x|) 


Q{\0t,T(v) -^D- 


Remark 3.2. Let us mention here the works of Knopova and Kulik IKK 14] . where the authors prove that the 
gradient of the density of a rotationally invariant stable process actually have the estimate: 


Pz(l,x) < 


1 


\d+a+l ' 


In other words, the exponent is actually bigger in this case. This fact allows to compensate the singularity of 
the gradient part when a < 1. This very important result is very specific to the case of the rotationally invariant 
stable process, as it is the particular form of the Levy measure that allows the growth of the exponent. In our 
general setting, we have been unable to prove a similar estimate, and had to restrict ourselves to remove the 
drift when a < 1. Observe nonetheless that this condition is relatively classic in the literature dealing with 
general stable process (see e.g. Kolokolstov |Kol00| 1. 

Consider now the integro-differential part of the kernel. For the small jumps part, once again, we observe 
that the operator acts on the variable x, and thus can be put on the density of the martingale. We use the 
representation in terms of symbols, denoting by <j> t (x,p) the symbol of an integro-differential operator $ 4 ( 2 ;, V^): 

®t(x, ^x)pm(t- t, (T{O t , T (y))- 1 (9t,T(y) -x)-^j 
( Z7r ) jR d 

Now, when $ t (x, V^) = L^(x, V x ) — L^(6 tt T(y), V x ), the small jump part of the difference of the generators, 
that is: 

(Lf(:r, V x ) - L“(0t, T (v), V*)) <p(x) 

(<p(X + z)~ ip(x) - (Vip(x), z)J l { | z |< (r _ t) i/c } (u(x, dz) - v(d t>T (y), dz)), 


I 

Jr 


m d 

denoting by l^(x, 0 — if 1 ( 9 t,T(y), C) the corresponding symbol, we have that: 

\C(x, 0 - if My), 01 < cs a My) - *r<“ Ai >icr. 
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Moreover, this quantity cf>t(x,—(a(6 t ,T(y)) 1 )*((T — t ) C) is smooth (in its </ argument) because of 

the truncation (see Theorem 3.7.13 in Jacob |Jac05| i. Consequently, 


T -t 

6 A 1 0 t , T (y) — a:| r? (“ A1 ) 


M x , 


-(^Ay^raT-tr^gT-tio 


is infinitely differentiable as a function of (j and uniformly bounded with all its derivatives. Therefore, it is in 
Schwartz’s space as well as its Fourier inverse. We have Vm > 1: 

(L™(x, V x ) - Lf (e t , T (y), S7 x ))p M (T - t, <r(e t , T {y))-\e t , T {y) - x) - $) 

^ A \ 9 t,r{y) - x\ ri{aM \ rr , |o-(6»t, T (?/)) _ 1 (^,T(2/) -x) -£| 

T-t [ ’ ^ + (T-t)Va 

Consequently, we recovered Lemma 2 in [SztlOl for the parametrix kernel, up to the additional multiplicative 
term (J A \0t.T{y) — x\ n '- aM ^) (T —f) _1 , which is the expected singularity for the kernel (see Kolokoltsov [KolOOj ). 
The upper bound follows from this upper bound and the control of the measure of the balls for P/v T _ t similarly 
to the derivation of the upper bound for the density, see Corollary 6 in [SztlOj and the proof of Theorem 1 
in ISztlOj . The upper bound for the small jumps part of the kernel follows. 

Finally, for large jumps, we see that the measure 1{|£| >{T-t) 1 / a }i v i x i d£) — v(O t ,T(y), d£)) is no more singular. 
Thus, we can write: 



J (p{t, T, x + y) - p{t, T , x, y)) l{|€|>(T-t) 1 /«}(^(a:, d£) - v(0t, T (y), <%)) 

< / p(t,T,x + £,y) - p(t,T,x,y) l{|£|>(T-tp/“} - H d t, T (y),d^)\ 

J R d 

< JA|x-0 t , T (2/)r ( “ A1) 


/S d_1 Jo 


+°° y( s ) 

p{t, T, x + sc, y)l{ s >(T-ty/°} dsyick) 


+ 7 f—p(t,T,x,y) ). 


For the last inequality, we exploited [H-5], We focus on the remaining integral term above. When the diagonal 
regime holds, the estimate is straightforward, as we can directly bound pit,T,x + s<;,y) < C(T — t)~ d ^ a < 
Cpit,T,x,y). The integral then yields the singularity (T — t) -1 . Therefore, we assume that 1 O t ,riy) — x\ > 
(T — t) 1 ' 01 . The regime of pit, T,x + sc, y) is given by |(9 tj r(y) — x — sc|- Thus, thanks to the triangle inequality, 
when | &t,Tiy) — x\ < 1/2s, or when s < 1/2| O t ^iy) — x\, the density p{t,T,x + sc, y) is off-diagonal with 
Pit, T,x + sc, y) < Cpit, x, y). 

Consequently, the problematic case is when s >c |0t,r(y) — x\. Indeed, in this case, pit, T, x + sc, y) can be in 
diagonal regime, whereas pit, T, x, y) is still in the off-diagonal regime. 

Assume first that [H-l -a] holds, and let us simply denote (c) the density of y on the sphere. Then, we 
have: 


r qis) dy 

/ / p{t,T,x + sc,y)l{ s >(T-t)V«}^T^ds— (c)dc 

J l/ 2 \ 6 t , T (v)-x\ Js*~ 1 s + ac 

p$/2\6t,T{.y)-x\ /• -/ \ j 

/ / p{t, T,x + SC, y)l { 5 >( T-i) v« } -p (c)s d_1 dsdc. 

Ji/2\e t , T (y)-x\ Js*-1 1 s +d dc 


Now, since s >c 1 9t,riy) — x|, we can take out of the integral. Also, the density ^(c) is bounded, so that 
we obtain: 
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/ 


3/2|e tlT (j/)-a:| 
l/2|f? t , T (y)-x| JS d 


[ P(t,T,x + s?,y)l {a > (T _ t) i/a } -^^ds^(<r)d<; 


c 


g(l^,Tfa) ~X|) 

|0t,r(2/) - x| a+d 7 0 /gd-i 


/*+oo 

JO 


p(t, T, x + SS, z)l {s >(T_ t) i, a} s d 1 dsdq. 


Finally, the remaining integral can be bounded by some constant as the integral of the density . Consequently, 
we obtained: 


M, T (y)-*\ r q(s) dp 

/ / P(t, T ,x +sq,y)l {s >, T _ t) i /a} - [ —ds—(<;)dq 

h/2\9 t , T {y)-x\ Js d ~i '« i+ “ <k 

T-t 


< c q(\°t,T(y) -x\) = c l 


I o t ,T{y) - x 


_ T \ot.+d 


T -1 1 d tiT (y) - x 


_ T \oi-\-d 


q(\°t,T(y) -x|), 


which is the off diagonal estimate for p when [H-la] holds, up to the singularity 1/(T — t). 

Now, assume that [H-l-b] holds. In this case, we can take out -Jry^ and integrate a density to get: 

r3/2\e t ,T(y)-x\ j- -f\ 

/ / p{t,T,x + sq,y)l {s> r T _ t) i/ a} -^j^dsp(dq) 

Jl/2\e t , T (y)-x\ Js d -i ^ 's 1+ “ 


< 


I i/2\ff t ,T(y)-x\ 

q(\0t,T(y) - x\) 


\0t,T(y) - x\ 1+a J 0 J S d -1 


r +oo 

J 0 


P y (t , T, x + s?, z)l {s > {T _ t) i/ a} dsn(d<;) 


< C 


q(\°t,T(y) - x\) 


\0t,r(y) - x\ 1+a ' 

Rewriting the right hand side to make the time dependencies appear : 
q(\0t,T(y) - x|) 1 ( T-t) 1+ ~^ 


\0t, T {y) - x\ x+a T-t\O t)T {y)-x\'+<* 


c 


1 (T-t) 


i+- 


T-t\e t , T (y)-x\^ 


In the last inequality, we recall that 7 < d, so that (T — t) ^ < 1. Now, we write: 

x 1 0 t , T (y) - x\ 1 ~ 1 q(\dt, T (y) - a;|). 


q(\0t,T(y) -®l) x (T-ty 


;<l(\0t,T(y) -x\). 


{T - t] ,+ " 


\0t,r(y) - z| 1+a ’ \9 t , T {y) - ^l Q+7 

Recalling we denoted by Q: 

Q(\0t,T(y ) ~ x|) = max(l, \e UT (y) ~ x|, 1 0 t , T (y) - x| 7_1 )g(|6> t ,T(2/) - x\), 
we finally obtain: 




T-t\6 t , T (y)-xy+° 


T -1 \0 t , T (y) - x\ a +*y 


Q(\0 t ,T(y ) -^D- 


In other words, we can correct the wrong decay by deteriorating the temperation. Consequently, the global 
upper bound for the kernel is the one announced. 
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To sum up, we deteriorated the tempering function in the following cases: 

• when the drift is bounded and [H-lb] holds. In this case, we replaced q by Q(s) = max(l,s 7 1 )q(s). 

• when the drift is unbounded and [H-la] holds. In this case, we replaced q by Q(s ) = max(l, s)q(s). 

• when the drift is unbounded and [H-lb] holds. In this case, we replaced q by Q(s) = max(l, s, s 7 1 )q(s). 
Note that when the drift is bounded and [H-la] holds, we do not need to deteriorate the tempering function. 

□ 


Remark 3.3. In the above proof, the temperation only serves to compensate the bad concentration in the 
generator. Also, we see that when the spectral measure p dominating the Levy measure v has a density on the 
sphere, then, the large jump part of the difference of the generators becomes: 



,T ,x + £, i/)l{|£|>( T _ t )i/ Q };/(d£) < C 



g(l£l) 

|£|<2+a 


<%■ 


Thus, when s | 0t,T(y) — x\, as in the last case discussed above, we have directly the good concentration index 

and the temperation is not needed. In particular, when q = 1, we recovered results in Kolokolstov IKolOOl . 

We have obtained the same type of estimate on the kernel and on the frozen density. Let us observe that the 
upper bound satisfies a ’’semi group” property in the following sense. 

Lemma 3.4. Fix r £ [f,T]. Let us denote 


Pc(t,T, x, y) 


(T - t)~ d / a 



|gt,T(y)-^| \ Q+ y 
(T-ty/* ) 


Q(C\O t , T {y) 


x|). 


Let Ci,C 2 > 0. For all r £ [t, T], there exists C 3 > 0: 



PC 1 ( t,T,x,z)pc 2 {T , T, 


z,y)dz < Cpc 3 (t,T,x,y). 


Proof. From the triangle inequality, we can write: 


\dt,T(y)-x\ < \0 t , T (y) - 0t,r{z)\ + \0t, T {z) - x\ 

< \0t, T (9 T , T {y)) - 9 t ,r{z))\ + 1 0t, T (z) - x\ 

< \0 r ,T(y) - z\ + \0t, T (z) -x\. 


The last inequality follows from the Lipschitz property of the flow. This means that we have either \0T,T(y) — z\ > 

(T-t)~ d '° 


C\0t, T {y) ^ x \ or 1 9t,r{z) — x\ > C\0t,T(y) — x\. Now, since for all t,T £ R+, u e-» 
decreasing, we have either: 


( 1 1 « ' 

\“ +7 

V (T-t) 1/a y 

I 


Q(Cu) is 


p Cl (t,T,x,z) < Cpc 3 (t,T,x,y ) or pc 2 (T,T,z,y) < Cpc 3 (t,T,x,y). 


Taking the corresponding density out of the integral, we can integrate the remaining density to one. The 
conclusion of the Lemma follows. 

□ 

We exhibit here some smoothing properties in time of the parametrix kernel. These properties will become 
crucial when investigating the convergence of the series (12.2311 on the one hand and the lower bound of Theorem 
ll.2l on the other. The following lemma is a regularizing effect in time of the parametrix kernel. 
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Lemma 3.5. There exists C > 1, co > 0 s.t. for all r £ (t,T), (x,y) £ (M d ) 2 : 

f 5 A\x-9 T ^{z)\^ aM) p{t,T,x,z)dz < C{r-t) u 

J R d 

[ 5 ^ T {y)-z\^ a ^v{r,T,z,y)dz < C(T - rf 

JR d 


As a corollary, we get that 


C dr [ \H(r : T,z,y)\dz<C[(T-tr+ (T-t) 1 ~ 1 / a l a>1 }. 

Jt J R d 

Thus, when integrated in time, the parametrix Kernel yields has a smoothing property in time. 

Proof. The two estimates are similar, we shall only prove one. Let us denote the quantity of interest: 

1= [ dzS A |S - 0 t , T (z) r (ttA1) . ( y ~T!'l <*+7 gfl* - M*)l)- 

jR d (i I *-0t,r(*) \ 

V + ) 


We split R d = D\ U D 2 , with 


D 1 = {2 £ K d ; \x - Ot, T {z)\ < C(t - t) 1/a } 
D 2 = {z£R d -\x-e t ^{z)\>C(T-t) 1 / a }. 


We write Ijo i for the integral over z £ Di. For z £ D i we have: 
(r - t)- d / a 


ft : l*-9t,T(*)h 
V + (r-t)V- ) 


a+7 


q(\x - 0 t ,r(z)\) <(T- t)~ d > a , \ X - 0 tiT (z) < (t - f)^ 1A1 /“). 


Also, Di is a compact and its Lebesgue measure is exactly (r — f) d /“, thus, we obtain I Dl < (r — t) , H 1A1 / a ). 
When z £ D 2 , we have: 


■y-d. 


Id 2 < 




< (r-t) 1 * 2 ? f 

J \z— 8 T ' t (x)\>C(T 


dz 


S A\z - e Tit (x)\^ aAls > 


_ t) l/c \z~0r,t(x) 1“+^ 

Observe that we used the Lipschitz property of the flow to switch from x — 9t,r{z) to z — 9 t ^{x). This allows 
us to change variables and set X = (z — 9 t j(x))/{t — t) 1 /®, we get: 


Id 2 <(j-t) 


i+2^ /' (r-t)^ 1 A «)|X|^“ A1 ) 

|a+7 


dX. 


J\x\>i \X\ C 

Thus, the result follows when a + 7 — d > r](a A 1). When it is not the case, we split again: 


6 A iz-e^tix)^^ 


dz = 


<SA|*-0 Tlt (a;)| , ' (aA1) 


\z— S T|t (x)|>(r—t) 1 /" \ z 9 T 't(x) | Q +2 J 1 >| 2 _ e rt (x)|>(r—i) 1 /" \ z 0T,t{ x ) | Q +'>' 

f £Alz~6> T ,t(z)| 7?(aA1) 


dz 


dz. 
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The second part of the right hand side is clearly a constant, bounding 5 /\\z — 9 T j{x)\ n ( aAl ^ < 5, since a + 7 > d. 
For the first part, we change variable again to Y = (z — 0 Tjt (x)), which yields when a + 7 — d < rj(a A 1): 

/ , . dY < C. 

J 1 > |Y|>(r—t) 1 /" l>T +7 

On the other hand, when a + 7 — d = T](a A 1) 

f 1 yn dY = [ lo g(l y l)](r-t)iA- < ~r\ lo g( r - *)l- 

Ji>\Y\>(T-ty/c \Y\ d ^ a 

Thus the proof is complete. 

□ 


3.3. Proof of Theorem ll.lt Uniqueness to the Martingale Problem 

We are now in position to prove the uniqueness to the martingale problem. Our approach relies on the 
smoothing properties of the parametrix kernel H. This method for proving the uniqueness to the martingale 
problem has initially been developed by Bass and Perkins [BP091 . and adapted to the parametrix setting by 
Menozzi |Menll| . For the sake of consistency, we choose to keep the notations in Menozzi [Menlll , as the proof 
is essentially the same. 

Proof. We focus on uniqueness. The existence can be derived from standard compactness arguments (see e.g. 
Chapter 6 in Stroock and Varadhan [SV79] , or Stroock [Str75| ). Suppose we are given two solutions P 1 and P 2 
of the martingale problem associated with L(-, V.), starting in x at time 0. We can assume w.l.o.g. that t <T, 
the fixed time horizon. Define for a bounded Borel function / : [0, T] x —> R, 


57 = E* 


f(s 1 X s )ds 


i e {1,2}, 


where (2T t ) t > 0 stands for the canonical process associated with (P*) ig { li2 }. Let us specify that S l f is a priori 
only a linear functional and not a function since P* does not need to come from a Markov process. We denote: 

S A f = S 1 f-S 2 f, 

and the aim of this section is to prove that S A f = 0 for / in a suitable class of test functions. Since this 
functional characterises the law of the process, we can conclude to the uniqueness. 

R), since (P*)ie{i, 2 } both solve the martingale problem, we have: 


if / e Cq , 2 ([o,t) x 


f(t, x) + r {ds + L s {x, Vi)/(s, X s )ds j = 0 ,16 (1,2}. 


(3.43) 


As a consequence we thus have that for all / £ Co’ 2 ([0,T] x 


S A ((d s + L s (xy a 


= 0 . 


(3.44) 


We now want to apply (13.4411 to a suitable function /. For a fixed point y £ and a given £ > 0, introduce 


for all / £ Cq’ 2 ([0,T) x R d ,R) the operator: 

V(t,x) £ [0,T) x R d ,G e,y f(t,x) = f ds f dzp s+e ’ y (t,s,x,z)f(s,z). 

Jt J R d 
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We define for all / £ Cq’ 2 {[0,T) x R d ,R): 

= [ ds [ dzL s {6 t , s+£ {y),V x )p s+e ’ y {t,s,x,z)f{s,z). 

Jt J R d 

We derive from the backward Kolmogorov equation for the frozen density that the following equality holds: 

d t G E ’ y f(t,x) + M E ; y f(t,x) = -f{t, x), V(t,x) £ [0 ,T) x R d . (3.45) 

Now, let h £ Cq’ 2 {[0,T) x R" d ,R) be an arbitrary function and define for all (t,x) £ [0,T) x 

<\> e ' v {t, x) := p t+e ' v (t, t + e,x, y)h(t , y), ^ e {t, x) := f dyG £ ’ y (0 £,y )(<, x). 

JR d 

Then, by semigroup property, we have: 

= f dy [ ds f dzp s+e ’ y (t,s,x,z)p s+e ’ v (s,s + £,z,y)h(s,y) 

J R d Jt J R d 

= [ dy [ dsp s+e ’ v (t,s + e,x,y)h(s,y). 

J R d Jt 

Hence, we can write: 

d t ^ e {t,x) + L t (x,S7 x )^ s {t,x) = j^ dy(d t G e,v (jf’ v {t, x) + Mtgcf> e ’ y (t,x)) 

+ f dy(L t {x^ x )G y ^y{t,x) -Mt;^it,x)) 


:= If{t,x) + IZ(t,x). 


Observe that from (13.451) . we have: 


If(t, x) = - p t+s ’ y (t, t + £,x, y)h(t , y)dy. 
Js. d 


Now, from Lemma m when e —>• 0 we have the convergence: 

f p t+E ’ y (t, t + £,x, y)h(t, y)dy —> h{t, x). 

J Rd E ->0 

Consequently, If(t,x) allows us to recover the test function h{t,x) when £ tends to zero, that is: 


lun\S A (I E 1 )\ = \S A h\. 


On the other hand, 


m^x) = f d dy(L t (x, V x )G v (t) e ’ v {t,x) - M E £ <j> s ’ y (t, xf) 

= J d d y J t ds(^L t (x, V x ) — L t {0t'S+e(y), V x )^p® +e,2/ (t, s + £,x,y)h(s,y) 

= dy dsH(t, s + £,x,y)h{s,y). 

Js. d Jt 
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From the controls of Subsection 13.21 specifically, Lemma [3~3l we have for all (t,x) £ [0, T] x R d : 

|7f(£, rc)| < |/i|oo /" dy [ ds\H(t, s + e, x, y)\ < C{T + e - t)“|/i|oo- 

J R d Jt 

Thus, denoting by ||S' A || := sup^i <x |S' A /|, we have: 


lim|S A (/f)| < ||5 A ||liminf|/I| 00 < C'||5 A ||(T-t)“|/ l | 00 . 


Now, from (I3.44|) with f(t,x) = \& e (t,x), we have 

5 a (( 5. + L(, V.))T e ) = 0 =► \S A (If)\ = |S A (Jf)|. 
Thus, for T — t small enough, 

\S A h\ = lim |S A (Jf)| = lim 15 A /|| < l/2||5 A |||/i| 00 . 

e->0 1 1 e->0 1 1 


By a monotone class argument, the previous inequality still holds for bounded Borel functions h compactly 
supported in [0, T) x R d . Taking the supremum over \h\oo < 1 leads to \\S A \\ < 1/2||S A ||. Since \\S A \\ <T — t, 
we deduce that |S lA | =0 which proves the result on [0,T]. Regular conditional probabilities allow to extend 
the result on R + , see e.g. Theorem 4, Chapter II, paragraph 7, in |Shi96] . □ 

3.4. Proof of Lemma 12.21 

In Subsection 13.II we have obtained estimates for both the frozen density and the parametrix kernel. In this 
section, we expose how these estimates are used to deduce the convergence of the parametrix series through the 
controls of Lemma 12.21 

Lemma 3.6. Fix t < r < T. There exists C > 1, tv > 0 such that for all (x, y) £ (M d ) 2 : 


p(t,r,x, z)H(t,T, z,y)dz < C (T - r) w + (r - f)“ + 


S A\6 t , T (y) -x\^ aA1) 


T-t 


p{t,T,x,y). 


Consequently, when integrated in time, we have: 

I P®H ( t , T,x,y)\ < c((T - t) u p(t, T, x, y) + p(t , T, x, y)^j , 

where we recall the notation p(t, T,x,y) = 5 A \0t,T(y) ~ x\ n ^ aA1 ^p(t, T, x , y). 

Proof. Let us recall that p and H are the upper bounds for p and H respectively, so that: 

\p® H(t,T,x,y)\ < [ dr f p(t,r,x, z)H(t,T, z,y)dz. 

Jt Js. d 


We investigate the integration in time. We have: 


p{t,T,x,z)H(r,T,z,y)dz = 


- P~ d /° 


(T-t) 


(1 i l*-0t,r(z)h 

V + (r-t)lTa ) 


a+7 


q(\x — 9t,r(z)\) 


S A\z-9r,T(y)\ v{aA1) (T-r)- d / a 

OL +7 


(1 , |z-0T,Tfe)h 

V + (r-r)i/- J 


T-t 


q(\z - d T , T (y)\)dz. 
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Assume first that \9t,r(y) — x| < C(T — t ) 1 / Q . The arguments differs according to the position of r in [t,T], 
First, assume that r G [2^,T], We have that r-txT-f, so that 


p(t,T,x,z) <(t -t) d/a x(T-f) d/a >zp(t,T,x,y). 

Consequently, we take p(r, x, z ) out of the integral and use the smoothing property of Lemma 13.51 

<5 A |.3-6»r,T(y)| r ' ( “ A1) (T-r)- d / a 


p{t,T,x,y) 


T-t 


a+7 


9(1* - 6> TiT (2/)|)dz 


/C , MrjMiV 
l 1 + (r-r)i7« J 

< C(T — r) aj_1 p(t, T, x, y). 


When, r € [0, we have T — r x T — t, and we have 


1 (T — (T — t\~ d / a l 

jT~“P( r i r, z, y) < T ] _ t < C—p(t, T, x, y). 


Next, we can bound 


(5 A |z — 0 T , r fo)| , ' (aA1) < Ct(6 A |0 t , T z - x^ 1 ) + S A |x - 0 t>T (y )| f?(Q!A1) ). 


Thus, we finally obtain: 


1 

T-t 


p(t,T,x,y) 


x(SA \0 t ,r{z) - 




( T -t) 


q(C\9 t , t(z) -x|) 


/C , |0 t , T (*)-*|V 

l 1 + (T-t) 1 '* ) 

x| 7(« a1 ) + s A |y tT ( y ) _ a:| , »( aA1 ))dz 


Assume now that |0t,r(y) — x\ > C(T — f) 1 /". In this case, the off-diagonal estimate holds for p(t,T,x,y), 
that is: 

p{t,T,x,y) X (r ^ -T^pr9(|0t,T(y) -x|). 

|0t,r(y) - x\ a+ T 

On the other hand, we have: 


1 0t,r(y) ~A < C T (\9tA z ) ~x\ + | 0 t ,t(v) -z|). 


In other words, we have either | 9 Tt r(y) — z \ > C\9 t ,T(y) — x|, or |0 tiT (z) — x| > C'|# ti T(y) — x|. Consequently, we 
split R d = T>i U T >2 with 


Di = {0 G M d , \9 T , T (y) - z\ < \9 t A z ) - z|}, 
D 2 = {z £R d ,\9 T , T (y) - z\ > \9 t , T (z) - x\}. 
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Now, when z G D\, we have that \6t,r{y) — a?| \9t,r{z) ~ x\, thus p(t,T,x,z) is off-diagonal and we can 

bound: 


p{t,T,X,z) 


< 


< 


(t - 

c ihhr^ mAz) 

I Ot,T(y) - x\ a+1 


x\) 

■ x \) ~P{t,T,x,y). 


For the last inequality, we used the fact that q is non increasing and that 7 + a > d so that the exponent in t — t 
is positive. Thus, we can take out p(t , r, x , z) of the integral, and use the smoothing property of H , Lemma 1531 
Denoting by Id t the integral in space in |p g) H\ where the space integration is over D 1 , we have: 


I„ < CHt,T,x,y) [ 

Jd, T-t (1 , \e T , T (y)-z\\ 

l 1+ (T-r)V- J 

< C'(T-rr- 1 p(t,T,x,y). 

When 0 G D 2 , we have |0 t ,t(2/) — z\ x |0t,T(y) — x|. In this case, observe that we have: 

g(| 0 t,T(y) - x \) 


1 -<r w ( T ^) 

~P\T — t, z, y) ~ 


-y(\ 9r,T{y) - z\) < (T — t)~ 


'' \e T Ay)-A a+1 \9t,r{y) ~ x \ a+1 

Thus, using Lemma [3~5l and recalling the constrain a + 7 > d, the integral becomes: 

r / ^ Q(\9t,T{v) - *1) ,rr ^ f j _ u , 

/Dl ^ C | 0 ,.r(i/)-»l^ (T - T ) • J^t, W ) 

x (<5 A |0 T , t (y) - ^“ A1 > + <5 A |0 t , r fo) - sc^C^ 1 )) 

9(l 0 t,T(y) - x|) 


C 


;(T-tY 


\0t, T (y) — x\ a+1 

x ((t-r) w +5A|^, T (i/)-ai|^ aA1 )) 

< c f( t - rr - + ^ l H ri T p( ° A1 > c.Tx,ri. 


T-t 


Thus the proof is complete. 


The following Lemma controls the second step of the iterated convolutions. 

Lemma 3.7. Fix t G (0, T\. There exists C > 1, u> > 0 such that for all (x, y) G (R d ) 2 : 

[ p(t,r, x, z)H(t, T, z, y)dz < c(jT - t) u ~ 1 + (r - T, x, y). 

J M 

Consequently, when integrated in time, we have: 

\p® H(t,T,x,y)\ < C{T — t) u p(t,T,x, y). 


□ 


Proof. The proof is similar to the previous one, but now, due to the presence of 5 A\d t , T (z) — a:| r? (“ A1 ) multiplying 
the first density, we do not use the triangle inequality anymore, because we are always in position to use 
Lemma 13.51 □ 
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3.5. Proof of the Lower Bound. 


Observe first, that due to the controls on the parametrix series, the convergence of the series actually yields 
a diagonal lower found for the density of (X t ) t> Q. Indeed, we have p{t , T, x, y) = pit , T,x,y) + p® H(t, T, x, y ). 
Also, we have the upper bound p(t , T, x, y) < p{t, T, x, y), which yields 

^ J /' ,SA\e T , T {y)-z\^ aA1) m w 

P® H{t,T,x,y) < ar p{t,T,x,z) -—- pir 1 T,z,y)dz 

Jt j R d i — T 

< (V - *r + A |0 t , T fo) - ^r (aA1) )p(t, T, x, y). 

Thus, in diagonal regime: \0t,T(y) — x| < C(T — t ) 1 / a , we have for f small enough p{t, T, x,y) > C(T — t)~ d N. 
In other words, we have a diagonal lower bound for the density of m- 

We now turn to the off-diagonal regime. The idea to derive a lower bound for the density is to say that in 
order to go from x to y in time T — t, we stay close to the transport of x by the deterministic system, for a 
certain amount of time, then, a big jump brings us to a neighborhood of the pull back of y by the deterministic 
system and the process stays in a neighborhood of this curve. 

In the off-diagonal regime: 1 0 t ^iy) — x| > C\T — i) 1 /®, we write from the Chapman-Kolmogorov equation 
for some to £ [t, T]: 


p(t,T,x,y) = / dzp(t,t 0 ,x,z)p(t 0 ,Tz,y)> / 

JVL d J B(6 t0 


,Tb/),C(T-t 0 )i/°) 


p(t,t 0 ,x, z)p(t 0 ,T, z,y)dz 


> 


> 


°(x to £B(O to , T (y),CiT-t 0 ) 1/a ) X t = x) inf p(t 0 ,T,z,y) 

\ / zGB(6t„ t(v),C(T— tn) 1 / 01 ) 


’(**0 € B(6 to , T (y), C(T - fo) 1/Q )|* t = x) C(T - t 0 ) 


zGB(e t0 , T (v),C(T-t 0 ) V“) 
-id\~ d / c 


Xt= x ]. To this end, 


Consequently we have to give a lower bound for P^X to £ B{9 t0t T{y),C{T — to) 1 ^) 

we introduce the process {X d ) t >o with jumps larger than <5 removed, for some <5 to be specified. Specifically, 
(Xf) s >o solves the SDE: 

X 5 s =x + J^ F(u,X s u )du + J a{u,X s u )dZt 

where (Z%) u>o is the process (Z u ) u >o with jumps larger that 8 removed. Its Levy measure is |^|<<5 

Now, observe that we can recover the process ( X s ) s >o from (X*) s > 0 by introducing the arrival times of the 
compound poisson process: 

x s = y] \z u i{\&z u \>8}- 

t<U<S 

Let us denote by (Tfc)fc>i the arrival times of the process (W)se[t,T] ■ We know that the variables Tk+i — are 
independent and have exponential distribution of parameter u(.B(0, <5)°). Then, we have: 

Vi < a < Tj, = X s s , 

X Tl = X* + <r(T u X* )AZ Tl , 

1 1 

VTi < s < T 2 , X s = X Tl + X s s - X*, 

J 1 

and so on. We refer to the Theorem 6.2.9 in Applebaum |App09| for a proof of this statement. We now split: 
p(x to eB(e to , T (y),C(T-to) 1/a )\x t = x') = P(X 6B(fl fc ,r(i/) 1 C(r-fo) 1 / |, );r 1 >io|^ = x) 


+r(x t0 e B(e t0tT (y),c(T- t 0 ) 1/a );Ti <to 


Xt = x). 
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Using the Markovian notations P*’ x (-) = P(-|X t = x), we thus focus on: 

P t ’ x (x to &B(e to<t (y),C(T-to) 1/a )-,T 1 <t 0 


= E* 


; (x 4o € B (0 to ,T (y ), C{T - to) l/a ) | T Tl ) l {Tl < to} 


where we denoted = er(Xf;s < Ti), the filtration generated by Xf until time T\. Now, by the strong 
Markov property, we have that 

P t ’ I (x (0 GB(^ 0 , r (y),C(T-f) 1 /“)|j ri ) = p T i> x Tj ^x to e B(9 to ,T(y),C(T — t 0 ) 1/a )) 


, B(8 t0 , T (y),C(T-t 0 y/°‘) 


p(T 1 ,to,X Tl ,z)dz. 


Thus: 


*'*(**„ e B{9 to , T {y),C(T < to) 


= E* ,x 




p(Ti,t 0 , X Tl , ^)dzl {Tl < to} 


Now, since X^ = X£_ + cr(Tj, X£,_ )AZt 1 , and since Ti is the first jump larger that <5, conditionally to X T - 
we have that <r(Ti, X T -)AZt 1 + X T - is a Poisson process on R d \B(0, S ). Thus, we have for all test function /, 
given X T -, the law of Xt x is: 

E[/(X Tl ) \X T - ] = E[/(a(Ti, X T -) A Z Tl + X T -)\X T -} 

= [ f(a(T 1: X )w + X ) ^ 

J{H><5} 1 1 P(S(0, 0 ) c ) 

Consequently, we obtain: 


E tx p(T 1 ,t 0 ,X Tl ,z)dz 

\JB(e to , T (y),C(T-to) 1 /°‘) 


X r, 


' B(9 tQi T(y),C(T—t 0 ) 1 /° l ) 


Tr J 

dz p(Ti,to,a(Ti,X^ r )w + X T ~ 5, z) 


'(dw ) 


i/(£(0,<S) c )' 


Now, we exploit the fact that Xj in independent and exponentially distributed with parameter v(B( 0, 5) c ) to 
write: 


' (X t0 e B(6 t0tT (y), C(T - Ti < to) 


= E 1 


rto 


ds 


dz 


'(dw ) 


L Jt 


B(e t0 ,T(y),c{T-t 0 ) 1 /o‘) Jwi d v(B(0, S) c ) 
p(s, to, *(s, X s s )w + Xf, z)u{B{ 0, S) c )e~ su ( B ^ S)C ) 


Observe that the quantity p(B( 0,<5) c ) gets cancelled. Now, we can give a lower bound by localizing the 
integral over w so that a(s, X d )w + Xf is close to 0 Stto (z). That is, where the density p(s, to, cr(X d )w + Xf, z) 
is in diagonal regime: 
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' (x t0 G B(e t0 , T (y), C(T - t 0 ) 1/a ): Ti < to) 


> E 4 ’* 


rto 


ds 


dz 


t JB(8 t0 , T (y),C(T-t 0 )«) J{\<r( s ,Xf)w+Xt-8 Stt0 (z)\<C(t 0 -s )«} 

xp(s, to, a(s, x s s )w + X s s , z) 


v(dw) 


> E 


rto . 

\ -d/a. f 

/ ds(t 

° - s) / 

.it ' 

> Jr 


dz 


xJb(<t(s, X s s )-\9 s , t0 (z) - X 5 S ),C (to - s) 1/a ) ) e-^ 0 ’ 5 ^ 


Additionally, we can lower bound the last probability by localizing X 5 S close to 9 Sjt (x): 

(x t0 G B(e t0 , T (y), C(T - t) 1/a ); Ti < to) 


> E* 


rto 


ds 1 


{|Af- 0 s , t (x)|<C( s -i) 1 /“} 


(to - s) 


—d/a 


dz 


’B(0t o , t (y), CiT-to) 1 /*) 

xv(B(a(s,Xl)-\6 sM (z)~Xl),C(to-s) 1/a ))e- s ^ B ^ . 


Now, from assumption [H-LB], is(B(0,5) c ) < 1 /8 a so that taking S = (T — i) 1 /" yields e S! '(- B (°>' 5 ) c ) > C. 
Also, since z G B(9t 0 ^(y), C(T — t 0 )“), by the Lipschitz property of the flow, 

O s ,t 0 (z) G 0 s , to (B(6» t0iT (y), C(T - to)*)) C B(O s , T (y), C(T - to)-), 

up to a modification of C for the last inclusion. On the other hand, Xf G B(0 Si t(x), s l ^ a ), thus, 

^XD-^M-Xl) G B(a(s,X 5 s )-\0 t ,T(y)- x), C(T -t) 1 ^). 


Now, from [H-LB], Ai ow is invariant under the action of a(t,x) for all (t,x) G R + x R. d , thus, if 


B(0 t , T {y)-x,C{T- 

we can use the lower bound in [H-LB] to get: 

v(B(a{8,Xi)-\9 8ito {z)-X s a ),C(to-s) 1/a )) > C (t 0 - s) 


(3.46) 


7/0 


Observe that this is exactly the condition (13.461) of Theorem 11.21 We thus obtain: 

Ft’^Xtc £ B(e t 0 tT (y),C(T -to) 1 /a y,T 1 <to) 

* '■y — d 

^ sl {IW 5 - e s,t(OI<C(s-i) 1 ''“} (to ~ s ) “ 

q(\a(s,X s s )-\9Uz)-X s s )\) - 

lB(8 t0 , T (y),C(T-t 0 y/°‘) Z \a(s,XS)-i(9 t , s (z)-XS)\'y+«_ 


> CE t,x 

X / 
JB 


rto 
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Consequently, since the function u i-a q(u)\u\ 7 Q is decreasing, the lower bound will follow from the upper 
bound: 


- X* a )\ < C\y - 0 T , t {x)\. 

We write from the ellipticity of a: 

\a(s,X s s )-\9 s , t0 (z)-X s s )\ < C\9 s , to (z) - X s s \ 

< C(\0 s , to (z)-6 s , t (x)\ + \6 s , t (x)-X s s \). 

Now, in the considered set, 1 9 s , t (x) — X %| < C(s — t ) 1/4 “ < C(T — t ) 1 /“ < C\9t,T(y) — x\. Thus, we have: 
\a(s,X 5 s )-\0 s , to (z)~X 5 s )\<C(\9 s , to (z)~9 s , t (x)\ + C\0 t , T (y)-x\). 

On the other hand, we can write: 


\0s,t o (z) ~ 8s,t( x )\ < \0s,t o (z)-O S ' T (y)\ + \0 StT (y)-Os,t(x)\. 

Thus, from the Lipschitz property of the flow, 

I Os,t(v) - 9 s<t (x )| < C T \O t ,T{y) - x\. 

On the other hand, we have 9 s ^{y) = 9 s ,t 0 ° 9t 0 ,T{y) so that: 

\°s,t 0 (z) - 9 8 , T (y)\ = \O s ,t 0 (z) - 9 S ,t 0 ° Ot 0 , T {y)\ < C T \z - 0 to , T (y)\, 

where to the get the last inequality, we once again relied on the Lipschitz property of the flow. We recall that 
|z — O to ,T(y)\ < C(T — f 0 ) 1/Q < | Ot,T(y) ~ x\, consequently we finally obtain: 


W(8,x s s )-\9 a , t o(z) - Xf)| < c T \e t , T {y) - X\. 

Using this last inequality to estimate the probability: 




> CE 4 


rt o 


C?Sl {| A 'f- e =,i(a:)l<(s-i) 1/ “} ( io “ S ) 


dz 


q(\9 t<T (y) - x\) 


' B(e t 0 , T {.y),c(T-t 0 y^ I e t,r(y) - z| 7+a 

> c(t - ^ 7 +I r ds P(|*f - 9 s , t (x)\ <(s- t ) 1 /“), 

\ 9 t,T{y) - a;| 7+ J t 


where (T — fo) d /“ comes from the volume of the ball B(9t 0 ,T(y),C(T — to) 1 /") obtained from the integral in 
dz. Using the diagonal lower estimates for the density, we actually see that P(|Xf — 9 s , t (x)\ < (s — t) 1/,a ) x 1, 
therefore: 


P 4 ’* (x t0 G B(9 t0 , T (y), C(T - t 0 ) 1/a ); Ti < to) 


>C(T-t 0 ) d /“(to-t) 1+2 ^ 


g(\ 0 t,T{y) - x\) 

\9t.r{y) — x \^ +a ' 
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Returning to the first estimate on the density yields: 


p(t, T, x , y ) > C(T - t 0 )- d ' a P*’ x (x t0 e B(O to>T (y),C(T - t 0 ) 1/a ); T x < t„) 

Z± g(\0t,T(y) ~ x \) 


> C(T-t 0 )- d/a (T-t 0 ) d/a (t 0 ~t) 1+: 


I OtAv) - z| 7+ ° 


= C(to-t ) 


l + z -d £(|6> t ,Tfa)-x|) 


| e t , T {y)-x^+ a ' 

Finally, to get the announced bound, we see that we have to choose to such that to — t x T — t. This gives: 

u+2 -d q(\0t,T(y) -a:|) 


p(t,T,x,y) > C(T-ty 


I O t ,T(y) - x\~ i+a ' 


which is the off-diagonal lower bound announced. 

Remark 3.4. We point out that the assumption [H-LB] appears quite naturally in this procedure as it serves 
here to give a lower bound on the u-measure of balls. Also, we see the that <j(t,x) preserves Ai ow is important 
so that we can use the lower bound for the u-measure of balls. 
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